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Abstract 

An integral in the sense of principal value of a singular function or of product 
of singular functions can appear itself as a singular function in some range of values 
of integration parameters. In this case, if necessary subsequently to integrate with 
respect to parameters, the problem arises about interpretation of the initial integral as 
a distribution over the integration parameters. A solution to this problem is offered, 
which is initiated by actual applications in quantum field theory. 

1 Introduction 

In many applications of quantum field theory one has to deal with integrals determined 
in improper sense. A classical example is the integrals of Feynman diagrams where the 
propagators of particles are determined (in Minkowski space) with a specific rule of the bypass 
of mass-shell singularities. With the loop corrections are switched on also the ultraviolet 
divergences appear in such integrals. The consistent and mathematically verified solution to 
the problem of "elimination" of the ultraviolet divergences is based on the conception of the 
extending of linear continuous functionals — what the coefficient functions of S'-matrix are 
— from a class of rapidly decreasing functions onto a class of arbitrary regular functions. 
As a matter of fact the solution is based on applying the distribution-theory methods El 
IU |H] . For the first time the idea of such an approach for solving the problem of ultraviolet 
divergences was put forward by N.N.Bogolyubov 6\ and afterwards was realized by him in 
collaboration with O.S.Parasiuk [J], and further by many other authors (see bibliography 
and the detailed account in monography [T]). 

A characteristic feature of the solution to the problem of ultraviolet divergences is the 
appearance of free finite parameters in the theory (in the renormalizable theories they are 
absorbed in favor of the renormalized constants). From the point of view of theory of 
distributions the mentioned property is entirely natural as it is connected with the operation 
of the extending of linear continuous functionals. Nevertheless, in many cases the solving by 
the same methods of other problems can be realized with the subsequent elimination of the 
ambiguities. Actually this occurs if the solution implies imposing of additional condition(s). 
Well-known example is the definition of a propagator in Minkowski space where the bypassing 
rule is fixed by the causality condition. In specific field-theory applications the elimination of 
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ambiguities can be carried out by the imposing of a self-consistency condition (see nontrivial 
example in [H]). 

Among the problems whose solution is based on the extending of linear continuous func- 
tionals but is realized without the emergence of ambiguities, a particular place is occupied 
by the problem of asymptotic expansion in a parameter of an integral defining some quantity 
(for example, the amplitude or probability of a physical process). Really, if the expansion 
is carried out before the calculation of the integral, then the expansion can led to the ap- 
pearance of singular functions in the integrand, which are not integrable in the conventional 
sense. In this case for giving a sense to divergent integrals one can take advantage of the 
theory of the extending of linear continuous functionals (the theory of distributions). How- 
ever, in contrast to the case of renormalizations, the ambiguities arising therewith must be 
completely eliminated — as the initial integral before performing the asymptotic expansion 
was well-determined (by our assumption) and its expansion contained no ambiguities. 1 A 
general recipe of the elimination of the ambiguities is explained in [D] . Nevertheless in some 
complicated cases that involve the calculation of the repeated integrals this recipe is insuffi- 
cient. For example, the method of |D| does not allow one to determine the expansion of an 
integral in a parameter if as a result the integral itself becomes a singular function in other 
parameters which, in turn, are considered as variables of the subsequent integration. The 
situation becomes even more complicated if the expanded integrand includes a product of 
singular functions. Exactly this occurs in the case of pair (multiple) production and decay 
of unstable particles when the process is described basing on the expansion in the coupling 
constant of Breit-Wigner factors which stand in probability (not in amplitude) jTTJ|. ITT] IT2*]. 
As was noted in jT2] i n this case already in the third order of the expansion the product of 
factors determined via the principal value (VP) emerges, and giving a sense to the product 
is a nontrivial problem. 

The purpose of the present paper is to investigate the above mentioned situation. In 
narrower sense our efforts will be directed on adding a sense to expressions of the type of an 
integral in VP sense of a singular function or of product of singular functions if the integral 
itself is a singular function of a parameter with respect to which a subsequent integration 
is supposed. The problem of elimination of the ambiguities arising at the determining of 
expressions of such type will be resolved by the imposing of a simple condition certainly 
satisfied in the actual applications standing beyond the present article. (The analysis and 
calculations directly in the framework of the mentioned applications see in |13j.) 

The structure of this paper is as follows. In the next section we find a solution to the 
problem of definition of integral in the case when the integrand contains a single singularity 
regularized by VP prescription but the integral itself is a singular function of integration 
parameters. Section |3] is devoted to the determining of an integral with the product of two 
VP in the integrand. In section0]we demonstrate the effectiveness of the designed method on 
a specific nontrivial example of calculation of integral. In section|5]we obtain a generalization 
of the method to the case of product of three, four and greater numbers of VP. In Conclusion 
the basic outcomes of the article are specified. 

1 The problem of asymptotic expansion under the symbol of integral arises when the integral cannot be 
calculated for technical reasons, but an expansion of the integral in a parameter is required. 
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2 Integral of a single VP-pole as a distribution over 
integration parameters 

Let us consider an integral of VP with varying (upper) limit of integration, 



y ^ 

F n (y)= [ dx VP— u(x, y) 
J x n 



Here u(x,y) is a weight function usually called a test function [2J El HI El- It is supposed 
that u(x,y) is distinct from zero only in a bounded area, and is finite and differentiable the 
necessary number of times. 

Through well-known formula for the principal value of a pole of degree n, 

x n (n-1)! dx n Vl 17 w 

supplied with the instruction to take an integral by the method of "integration by parts" 
[HI IU E] j one can obtain 



F n (y) ' 



n-1] 



J fc =i y 



(3) 

Here u^ k '(x, y) = d k /dx k u(x, y) and u^ k '(y, y) = d k /dx k u(x, y)\ x = y - 

Expression (JBJ) defines function F n (y) at any y ^ 0. At y = the value of F n (?/) is not 
determined as at n = 1 expression contains a logarithmic singularity, and at n > 1 a 
power one. Correspondingly, function F n (y) at n = 1 is integrable in a neighborhood of 
y — 0, and at n > 2 is not such in the conventional sense. Nevertheless, the meaning of 
function F n (y) can be extended to make it integrable in the sense of distributions. The 
necessity in this operation arises if F n (y) should be further integrated with respect to y. 

As a recipe, a solution to the problem of extending the meaning of a function containing a 
pole consists in assigning a principal value to this pole and adding a functional concentrated 
at the point of singularity of the pole. Actually the mentioned functional must be a sum of 
Dirac delta function and of its derivatives with arbitrary coefficients, and the degree of the 
high-order derivative should be less by 1 than the degree of singularity of the pole [3| (the 
degree of the high-order derivative is fixed by the number of necessary subtractions in the test 
function after which the integral becomes well-determined in the conventional sense). In the 
case of formula Q the solution consists in replacing each l/y k by VP (l/y k ) + J2 k =o Q 5^ l \y), 
where Ci are the coefficients describing the parametric ambiguity. 

Now we show that the ambiguities in formula (JHJ), arising at the extending of meaning 
of the poles, can be completely removed by the imposing of a condition of independence 
of the result of integration from the order of calculation of repeated integrals. In effect, 
the mentioned condition means equivalence of the result of the repeated integration to the 
result of the multiple integration. (Of course, this is not the only possible condition, but it 
naturally arises at solving enough wide class of problems.) 
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So, let us turn to the initial formula (JTJ) and consider functional F n [u] specified by repeated 
integration, 



dy 



y 1 

da: VP — u(x, y) . 

^ ry.1l V ' ^ ' 

-OO X 



(4) 



Our aim is to determine this functional so that to make it be equal to a multiple integral, 



da; dy VP — 8(y — x) u(x, y) 



x' 



(5) 



or to a functional determined in a similar way, but with opposite order of the repeated 
integration, 



r 

Fn[u] = / 



oo ] 

dx VP 

oo X 



it 



dy u(x,y) 



(6) 



It should be emphasized that in contrast to (j3J) functional © is well-determined since the 
integral dy in © defines a function from the space of test functions of one variable if u(x, y) 
belongs to the space of test functions of two variable. 

By making use of the definition (0) for principal value and carrying out obvious calcula- 
tions, we obtain 



F n [u] 



n-l) 



dx \n(\x 




n—l 



u(x,y) - 



fc=0 



dx k 



d 



n—k—l 



dx 



n—h 



- u(x,y) 



(7) 



In all integrals are considered in a proper sense and converge. Therefore we can change the 
order of integration in the first term. In the second term we substitute y for x. Then we get 



T 

J n 




dx ln(|x|) 



dx T 



n-l 



-U[X 



,y)-Mly|)ETT 



d< 



n—k—l 



dx 



n—k- 



- u(x,y) 



\x=V. 



With taking into consideration (0) the later formula may be rewritten as 



^n[<j>] 



1 



(n-l)\ 



dy { j dx ln(\x\)—4>(x,y) 



—oo 

nn-1 n-l i Q 

- H\y\)^— 1 ^y)\,=y + ^ k - i y-y p 

k=l 



n—k—l 



yk Q x n-k- 



{x,y)\x= 



(9) 



2 The introducing of formula JSJ may be considered as a heuristic trick useful for the transition to formula 
©• Nevertheless, a precise mathematical sense can be added to as well. Really, since a singular 
distribution can be considered as an improper limit of a conventional function with respect to some parameter 
(by means of auxiliary regularization) , the multiple integral in JSJ can be considered as a conventional integral 
with the posterior transition to the limit. Formulas and 10 are equivalent from this point of view. At 
the same time, functional Q remains uncertain as long as the transition to the limit (removal of auxiliary 
regularization) before the calculation of integral dy leads to a meaningless result. In effect wc determine 
functional by the imposing of condition to remove auxiliary regularization after the calculation of all 
repeated integrals. 
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It can easily be seen that the expression in the curly brackets in © coincides with that in 
the r.h.s. of formula ©, considered with substitution VP-poles for the ordinary poles and 
without the adding of ^-functions and its derivatives. 

So, we have shown that the condition of independence of the result of the repeated 
integration of VP from the order of calculation of the repeated integrals (the requirement 
of equivalence between repeated and multiple integrations) implies the VP prescription for 
the poles arising at carrying out the first integration in the repeated integrals. In the case 
when both limits of integration are variables and when VP is appeared with the "shifted" 
argument, one can easily obtain the following general formula: 



dx VP- 



-y) n 



u(x, 



[n- 



dx \n{\x-y\) w (n) (x, ...)+ ln(|6-y|) u^\b, . . .) - ln(|a-j/|) u^a, . . .) 



n-l 



-EC*- 1 )' 



k=l 



VP 



(b-y)> 



■ u 



(ra-fc-l) 



- VP 



(a-y)> 



, u {n-k-l) f 



(10) 



Here u(x, . . .) = u(x, y, a, b) and the superscripts in brackets mean the partial derivatives of 
the corresponding order with respect to the first argument. In the case when an integration 
dy is supposed as the next one, it is convenient to represent formula (fTUJ) in the form 



dx VP- 



[x- 



n-1 



-y) 1 

■ fc=0 



U[X 



d k 



n-l)! 

a 



dx ln(|x— y\) w (n) (x, . . .) 



■V\)u 



dy' 



\n(\a-y\)u 



(n-fe-l)/ 



Here the derivatives with respect to y are understood in the sense of distributions, i.e. they 
are to be moved by the rule of "integration by parts" at the next integration dy. 

Formulas (fTU|) and (fTT|) admit various writings in special cases. For example, if the 
dependence on y in the test function can be separated out into a factor (j)(y), then there is 
the following formula: 



dx VP- 



x — y) n 



u(x,a,b) <j>{y) 



n-l 



dy n 



dx ln(|x— y\) u(x, a, b) . (12) 



Here the derivatives are again understood in the sense of distributions. The validity of (|12p 
follows from the comparison of the result of applying (fTTj) to the l.h.s. of (|12|) with the result 
appearing in the r.h.s. after the change of integration variable x — > x + y and carrying out 
then the direct calculation as in conventional integral. 

The other important corollary of formula (fTU|) appears in the case when the test function 
does not depend on variable x within the limits of integration (a ... b). In this case at n = 1 
we get 

ft 

dx VP u(y, a, b) = u(y, a, b) 



x 



y 



\n(\b-y\)-\n(\a-y\) 



(13) 
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At n > 2 we obtain 



dx VP 



(x - y) n 



u(y,a,b) 



n — 1 



VP 



ip-y) 



n-l 



VP 



(a - y)^ 



(14) 



Let us remember once again that in formulas (fTU|) - (fT^ the quantities a, b, y (or, at least, 
some of them) are considered as variables of a subsequent integration but not as parameters. 
In the case when the mentioned quantities are considered as parameters, the symbol VP in 
the r.h.s of the above formulas should be omitted and the formulas themselves have a sense 
only at uncoincident arguments occurring in the logarithms and/or poles. 



3 Product of two VP 

Now let us consider an integral with a more complicated structure containing a product of 
two VP-poles: 

f 1 1 

/ dx VP- — VP- — u(x, ...). (15) 

a 

Unfortunately, the formulas of the previous Section do not allow one to add a sense to 
expression (|TK|) since at z\ = z 2 the integrand containing a product of singular functions is 
undetermined. Nevertheless, the method of the previous Section, basically, can be applied 
in this well. Really, we can at first determine integral (fTB^) on the assumption that 

z\ and z 2 are the parameters of integration, not equal one another. As a result we obtain 
a singular function of Z\ and z 2 . Then we extend the meaning of this function in the sense 
of distributions. For removal of ambiguities originating therewith we proceed to a repeated 
integral of triple multiplicity containing integration with respect to z\ and z 2 , and determine 
this integral by the imposing of condition of independence of the result from the order of 
integration. (It should be used the property that the integral is well determined if at first 
the integration is carried out with respect to z\ and z 2 , and only then with respect to a;.) 

However, the above mentioned method in the case of two poles is found too cumbersome 
and hardly probable is justified, especially if there is a greater number of VP in the integrand. 
So let us use of a trick based on an independent determination of the product of two VP 
poles. (We emphasize, once again, that the necessity of reference to integral (|T5J). instead 
of at once to the triple repeated integral with another order of integration, is caused by the 
fact that the indicated order of integration provides a practical solvability of a problem in 
some applications.) 

At first we consider the case n\ — n 2 — 1, and at x ^ z±, x ^ z 2 , Z\ ^ z 2 examine the 
following formula with the ordinary poles realized in the sense of conventional functions: 

1 11 

X — Zi X — Z 2 Z\ — z 2 

The expressions in the both sides of formula (|T6|) may be also considered as functionals 
determined on space of test functions vanishing at the coincidence of any pair of arguments. 
On such space of test functions there is still an equality ()16|) between both functionals. 



X — Z\ X — Zo 



(16) 
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Now let us state a problem about the extending of these functionals onto the all space 
of test functions. The solution we realize in two steps. At first we determine each pole 
by assigning VP and adding the ^-function with arbitrary factor. Then we determine the 
product of VP-poles. For solving the latter problem we use the remarkable property of 
formula (j!6|) which consists in the fact that its l.h.s. may be determined as a product of two 
VP if at first the integration is meant with respect to z\ and/or z 2 and only then with respect 
to x. At the same time, the r.h.s. of (fTBj) with VP attached to the poles is well determined 
if, conversely, the integration with respect to x at first is supposed. Let us call the order 
of integration by regular if it coincides with the above mentioned one, and by irregular 
otherwise. Then we see that the r.h.s. of formula (fTp^l can be used for determining the 
l.h.s. in the case of irregular order of integration, and the l.h.s. can be used for determining 
the r.h.s. in the corresponding case. In both cases the equating of one side of the formula 
to another side should be treated as the extending of a linear continuous functional. As 
was directed above, this operation is not unambiguous and can make sense only up to a 
functional concentrated at the point of singularity (uncertainty) of the initial unextended 
functional. In our case this means the necessity of adding the product of two 5-functions 
with arbitrary coefficient to one of the sides of the resulting relation. 3 Below we present the 
result obtained with the taking into consideration of the symmetry at reading the formula 
from left to right: 



VP- 



VP- 



x 



X — Z 2 de f 



VP- 



1 



zi - z 2 



VP- 



1 



VP- 



x 



Zl 



Z2 



(17) 



-VP- 



Ci 



S(x — zi) — 5(x — z 2 ) + C 2 S(x — zi)5(x — z 2 ) . 

Zi - Z 2 L 

Here the following notes are in order. First, we have not shown the contribution of 5- 
function that has to be added to VP {z\— z 2 )~ x in the r.h.s. since this contribution is zero in 
view of nulling the expression in square brackets at z\ = z 2 . Second, we have not written the 
contributions of 5-functions added to single poles in the l.h.s. since by virtue of the symmetry 
the structure of the relevant contributions coincides with the structure of the second term in 
the r.h.s. of the formula. (Therefore, the mentioned contributions are absorbed by the second 
term in the r.h.s.) Third, we consider the product of (^-functions in the r.h.s. of (j!7j) to be 
determined so that the result of its integration is independent from the order of calculation 
of integrals. The latter requirement implies that 5(x — zi)5(x — z 2 ) = 5(zi — z 2 )5(x — z 2 ) = 
5(zi — z 2 )5(x — zi). Under this condition formula (fTTj) can be read in the opposite direction, 
from right to left. In the latter case all terms with (^-functions are to be transferred to the 
l.h.s. of the relation. 

Coefficient C\ in (|17|) can be determined through condition of the recursive invariance of 
the formula. Namely, let us demand the invariance of the formula at the removal of square 

3 A general theory of extending of linear continuous functionals (regularization of singular functions) 
is discussed with enclosing of numerous illustrations in [3] and In the context of the problem under 
consideration it is appropriate to carry out an analogy with the theory of renormalizations in the field theory: 
the replacing of the ordinary poles by VT'-poles and the adding to them of 5-functions can be compared 
with the elimination of UV divergences accompanied by emergence of the ambiguities in subdiagrams; the 
determining of the product of two VP and the adding of the product of <5-functions corresponds to the 
elimination of UV divergences and emergence of the ambiguities of the overall type in the diagram (see 
Section 29 in monography P). 
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brackets in the r.h.s. and determine the arising products of VP by means of the formula (fTTj) 
itself. Then we get 

Ci = 0. (18) 

Unfortunately, in doing so the coefficient C 2 remains uncertain, i.e. it cannot be fixed by 
the reasons of symmetry. So, the determining of C 2 is only possible by the imposing of an 
additional condition. As such condition we consider the requirement of independence of the 
result of integration from the order of calculation of the repeated integrals of the r.h.s. of 
()17|). In the case of irregular order of integration this requirement will fix C 2 by equating the 
result of integration to that obtained at the regular order of the calculation of integrals. The 
easiest way to actualize this condition is to carry out calculations at some special choice of 
the test function. The elementary choice is the unit function in the integral with finite limits. 
So, let us consider the following repeated integral of triple multiplicity: 



111 



1, 



[dx I dz a [ dz 2 VP—*— VP— *— 

J J J X — Zi x — z 2 



(19) 



000 

Its direct calculation with the use of (|T3"j) leads to result 

On the other hand, again by direct calculating we can get 
111 



(20) 



1zi,z 2 ,x = J ' dzi J dz 2 J dx VP 
000 

From (j201), (j2U) and (H?|) we conclude 



1 



zi ~ z 2 



VP- 



1 



-VP 

X — Z\ x — z 2 



7C 2 . 



(21) 



C 2 = n 2 . 



(22) 



Result (|22j) can be obtained also on the basis of some skilful manipulation with Sokhotsky 
formula. Really, let us consider a product of two simple poles with the "causal" bypass of 
singularity (see below, formula (|23))). It should be noted that this product is well determined 
on the space of test functions under consideration [Sj. If integration dzi and dz 2 at first 
is implied (with any test function), then the Sokhotsky formula can be applied to each 
multiplier. Moreover, in the resulting expression the brackets can be removed, as well. 
Eventually we obtain 



1 



1 



x — z\ + 1O x — z 2 + 1O 
1 



1 7T 



VP- 



VP VP- 

X — Z\ x — z 2 

1 



X — Z\ 



5(x - z 2 ) + VP 



x - z 2 



5(x — zi) 



7T 2 8{x - z x )8{x - z 2 ) ■ (23) 



Here the point placed above the equality symbol recalls that the equality has a sense only 
at the particular order of calculating the repeated integrals. 
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On the other hand, provided that the integration dx at first is carried out, the expression 
in the l.h.s. of (J2*3)l can be transformed to the form (see the first note below formula Q17])) 



VP- 



1 



1 

x — Z\ 



1 

z 2 + i0 



Z\ — Z 2 IX — Z\ + lO x 

Again by applying Sokhotsky formula and removing the brackets we get 
1 1 



(24) 



x 



Z\ + lO X — Z2 + lO 



VP- 



1 



- ^2 



VP- 



- VP- 



1 



X — Z\ 



x - z 2 



mVP- 



Z\ - z 2 



S(x - z x ) - S(x - z 2 ) . (25) 



The point under the equality symbol in ()25|) recalls that here another (particular) order of 
calculation of the repeated integral is supposed. 

Now by virtue of independence from the order of calculation of the repeated integrals of 
the product of "causal" factors, we equate the r.h.s. in (|2Hjl and (|25jl . Then again we obtain 



So, ultimately formula (jl7j) takes the form of 



VP- 



VP- 



VP- 



x 



Z 2 



Zl - z 2 



VP- 



X — Zi 



VP- 



x - z 2 



+ ix 2 5(x - zx)S(x - z 2 ) . (26) 



The meaning of this formula is as follows: its r.h.s. defines the l.h.s. in the case of irregular 
order of calculation of integrals, at first dx and only then dzi and/or dz 2 . This definition 
provides the equality of the result of integration of the r.h.s. to that obtained at the regular 
order of integration of the l.h.s., at first d;?i and/or dz 2 and only then dx. (It should be 
emphasized, once again, that the imposing of another condition may change the value of C 2 .) 

A generalization of formula (|26j) to the case of the product of poles of an arbitrary degree 
can be easily generated through the use of relation 



VP- 



1 



1 712-1 



VP- 



{x-zx)^ [x-z 2 y* (th-1)!(ti2-1)! dz? 1-1 dz\ 
By substituting (|26J) into ()27|) we obtain after simple calculation 
1 „„ 1 



VP- 



VP- 



x — Zi 



x - Z 2 



(27) 



VP 



(x - Zi) n ^ 



VP 



7T 



(x - Z 2 ) n * 
ni~ 1 

+ e ( n2+ r 



k=0 
U2 — 1 



( ni -l)!(n 2 -l)! 

1 



VP 



{zi 



+ E ( ni+ fe fc_1 ) (-) fe vp 



z 2 

1 



\n 2 +k 



VP- 



1 



712-1) 



(x - z 2 ) 



[X 



k=0 



(z 2 ~ Zi) 



rii+fc 



VP 



Zl 

1 



mi— k 



(x - Z 2 ) n2 ~ k 



(28) 



Returning to integral (J15)) introduced in the beginning of this Section, we see that it 
can be determined by means of formula (J28j) and then can be calculated with the aid of the 
formulas of the previous Section. 
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In conclusion of the present Section we note that formula ()26|) in effect is not completely 
new. In particular, in ^1] a similar formula was derived, presented in the form of a relation 
between repeated integrals of double multiplicity, and in [TS] an equivalent up to notation 
formula was presented (but without derivation and comments). However the derivation 
in |14j . which is based on calculation of the conditional limits with respect to parameter 
of the conventional integrals, substantially differs from our derivation which is based on 
the extending of linear continuous functionals. The basic advantage of our derivation is the 
universality and flexibility of the mathematical tools in operation. This reflects, in particular, 
in a possibility of automatic generalization of the results to the case of any multiplicity of 
the integrals, and also in the extreme transparence and brevity of the proposed solution. 



4 Nontrivial example of calculation of integrals 

Let us consider an example of the use of formula ()26|) close to that which appears in some 
actual applications. Namely, let us consider an integral over a simplex 

°° 11 
I( z ) = Ax Ay 9(2 + z-x-y) VP—- VP . (29) 

y 

At once we note that at z > the point of singularity of both poles {x — 1, y = 1} certainly 
falls on the integration area. At z < only one of the poles can be singular. The case z = 
in some sense is transitional. Simultaneously this case is specific-singular because at z = 
the ambiguity of ^-function is joined to the singularity in integrand. (So, at z = integral 
(|2HJ) requires of additional determining. A possible way is noted in the footnote on page (jlj).) 

Integrals of the type of (J2"9*j) arise at calculating the probabilities of the processes of pair 
production and decay of unstable particles in the approach of a modified perturbation theory 
(MPT) based on expansion in the coupling constant of the Breit-Wigner factors standing 
in the probability (not in amplitude). Both V^P-poles in (J29|) in this connection correspond 
to the particular contribution emerging from the product of two Breit-Wigner factors in 
the next-next-to-leading order of the expansion of the cross-section. The points x = 1 and 
y = 1 correspond to the positions of the mass-shells of resonances. Quantity z stands for the 
energy of exclusive process counted off from the threshold of the pair production. (See [TJ] 
for the establishing of correspondence, and 13J for computation of an actual processes.) 

Let us represent (J2lj|) in the form of repeated integral and take advantage of formula (fT5|) . 
Then we obtain 

2+z 2+z-x 2+z 

I(z)= [AxVP / Ay VP = / Ax VP ln|l+z-d. (30) 

J x — I J y — 1 J x—1 



Unfortunately, at z = the last integral in (f3T)|) is not determined. Nevertheless at z ^ 
its calculation can be carried out in a direct way by separating the range of integration onto 
the sub-ranges. (Note that in some complicated cases this is not always possible to make.) 
Omitting tiresome calculations, we write down the result at z > 0: 

2 

I(z) = 2 dilog(l + z~ l ) + ln 2 (z) - ^ , (31) 
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dilog(z) =f'dt^. (32) 

The more perfect method of calculation of I(z) is based on the change of variables 
X + U = £,, x ~y = 2v (therewith the symmetry of the going through of the integration area 
is achieved) and on the usage of formula 

2+2 ?/2 1 

d£ / dr? VP — VP- 



-f/2 
2+z t/2 



77 + <e/2 - 1 77 - e/2 + 1 



/ d£ y dr/JvP 



- 5 /2 
2+2 



^-2 



VP ; VP 



r/ + e/2-l 77-^/2 + 1 



2 / de ln 5 ^ -n 2 6(z). (33) 



The last integral in (}3*3*j) can be calculated at any z. In particular, at z > — 1 we get 

7T 2 

I{z) = -2 dilog(l + z) + — - 7T 2 6(z) . (34) 

At 2 > the expressions in the r.h.s. of (}3~T|) and (J3*lj) are equal each other by virtue of 
the relation 

2 

2 dilog(l + z~ l ) + 2 dilog(l + z) + ln 2 (z) + y = . (35) 

The validity of ()35|) can easily be verified by the differentiation of the l.h.s. with taking into 
consideration (|32|) and by calculating separately its particular value, for example, at z = 1. 

Thus, at z > the both above mentioned calculations are equivalent. Nevertheless, the 
second method of calculation of integral (|29|) allows one to solve the problem of the going 
through the range of the threshold of pair production of unstable particles, a stumbling- 
block from the point of view of calculations of jT^j. (In actual applications a singularity 
arises at z — 0, which needs in a regularization.) Besides, outside the threshold (at z ^ 0) 
the calculation by the second method is considerably simpler, which is very important from 
the point of view of the practical solvability of a problem. 



5 Product of several VP 



The results of Section El may be generalized by induction to the case of any number of VP- 
poles. So, the product of three poles of degree 1 can be determined via multiplying both sides 
of formula (j2*U|) by one more pole. As a result in the r.h.s. we obtain the product of no more 
than two VP-poles of x. By virtue of ()26|) this product is well determined. Furthermore, by 
sequential applying the result can be reduced to the form of a sum of single poles of x. 
Written in the completely symmetric form, the result is 



nvp 



71=1 



3 

E 

71=1 



VP- 



nvp 



7T~ 



Y[ 5(z n - Z k ) + 7T 2 



Y[ 6(x - z k ) 



(36) 
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Formula for the product of four poles can be obtained by multiplying both parts of (|36p 
by one more pole, or by multiplying formula (|2*Bj) by itself. Both methods after reducing to 
the form of symmetrized sum of single poles of x lead to result: 



4 1 4 1 

n=l x ^ n n=l x / " r ' 



Here indices and / run values 1,2,3,4 except for the values indicated under the symbols of 
sum or product. The summation in the second term in the r.h.s. is carried out over the all 
permutations of {zi, z 2 , Z3, Z4}. 

The process of multiplying by a simple VP-pole can be continued. In doing so at each 
step of the induction in the r.h.s. we obtain the product of only two VP-poles of x, which 
in view of (|26|) is a well-determined quantity. After the necessary number of steps we can 
obtain a formula that determine the product of any number of simple VP-poles. Then by 
the differentiating with respect to parameters, similarly as in ()27|) . we can derive the result 
for the product of any number of VP-poles of any degree. In effect this procedure is trivial. 
In view of awkwardness we do not write down the ultimate result. 



6 Conclusion 

So, if an integral in the sense of principal value of a singular function or of product of 
singular functions appears by singular function of the parameters of integration, then the 
integral can be determined in the sense of distributions. The ambiguities arising therewith 
are in control and can be eliminated by the imposing of the condition of independence of the 
result of integration from the order of calculation of the repeated integrals. (The mentioned 
condition, of course, is not unique but naturally arises at solving the wide class of problems.) 
In the case of a single pole with VP-prescription in the integrand, the result is described by 
formulas (J1(J|) and (|11|). The essence of these formulas consists in assigning VP-prescription, 
again, for the poles arising at the conventional calculating of integral. The case of product 
of two VP-poles in the integrand is nontrivial, but it can be reduced to the case with a single 
VP-pole via the reduction formulas (|2*Bj) and (}2*H|) . The case with many VP-poles is easily 
considered by induction. 

The results of the present article are extremely important for a systematic description 
of the processes of pair (multiple) production and decays of unstable particles in the high 
orders of perturbation theory. Furthermore, in view of the obvious universality they can 
be applied in other applications, as well, containing the repeated integration of singular 
functions determined in the sense of principal value. 

The author is grateful to V.A.Petrov for valuable notes, and also to A.I.Alekseev for the 
indication to Ref.|15j. and A.Bassetto for the indication to Ref . |14j . 
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